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I.  INTRODUCTION 

Observational  studies  of  the  Gulf  Stream  front  have  shown  the  horizontal 
separation  between  the  Gulf  Stream  surface  and  subsurface  (15°  isotherm  at 
200m)  fronts  to  be  quite  variable.  Early  studies  such  as  those  by  Hansen  and 
Maul  (1975)  and  Khedouri,  et  al .  (1976)  had  insufficient  wind  data  to 
definitely  correlate  changes  in  separation  with  the  passage  of  storms, 
although  Hansen  and  Maul  did  suggest  that  there  was  a  relationship. 

During  August  1981  the  Naval  Oceanographic  Office  (NAVOCEANO)  conducted  a 
survey  of  the  Gulf  Stream  northeast  of  Cape  Hatteras  to  study  the  relationship 
between  the  paths  of  the  Gulf  Stream's  surface  and  subsurface  fronts  and  the 
local  wind  stress  (Horton,  1982).  During  this  experiment  the  center  of 
Hurricane  Dennis  passed  from  west  to  east  approximately  30km  to  100km  south  of 
the  portion  of  the  Gulf  Stream  under  observation.  At  the  same  time  the 
average  separation  between  the  surface  and  subsurface  fronts  increased  from 
about  1km  to  about  19km.  Since  the  zonal  component  of  the  winds  was  from  the 
east  during  the  passage  of  Dennis,  the  increase  in  separation  is  consistent 
with  a  displacement  of  the  surface  front  to  the  right  of  the  wind  by  advection 
in  the  mixed  layer. 

In  this  study  we  model  and  predict  the  displacement  of  the  surface  front 
of  the  Gulf  Stream  due  to  short  but  intense  wind  events  and  compare  the 
results  of  this  effort  to  the  observed  displacement  of  the  surface  front  by 
Hurricane  Dennis.  The  model  is  a  relatively  simple  extension  of  a  one¬ 
dimensional  bulk  mixed-layer  model  to  two  dimensions  to  permit  the 
representation  of  a  cross  section  of  the  mixed  layer  in  the  frontal  zone  of  a 
strong  current.  No  downstream  variations  of  the  front  are  allowed.  A 
secondary  purpose  is  to  illustrate  how  strong  mixing  or  cross  stream  mass 
convergences  might  lead  to  the  formation  of  new  surface  fronts  or  the 
destruction  of  existing  ones. 

A.  Model  Background . 

Earlier  efforts  by  de  Szoeke  (1980)  and  Cushman-Roisin  (1981)  have 
extended  one  dimensional  mixed  layer  models  to  two  dimensions,  but  their 
efforts  differ  substantially  from  this  one  in  that  their  initial  states  do  not 
contain  fronts.  These  models  start  with  meridional  temperature  gradients 
characteristic  of  the  open  ocean.  They  are  primarily  interested  in  the 
formation  of  horizontal  variation  (especially  fronts)  in  the  surface 
temperature  field  on  the  seasonal  time  scale,  although  the  de  Szoeke  model  is 
also  valid  during  the  first  inertial  day. 

According  to  Cushman-Roisin  (1981)  there  are  two  major  classes  of  mixed 
layer  models:  a)  turbulent  erosion  models  (TEM),  and  b)  dynamic  instability 
models  (DIM).  In  order  to  describe  the  differences  between  them,  we  write  the 
simplified  potential  energy  (PE)  budget  for  the  mixed  layer  as 

4^  =  F  -  D  +  E  .  (!) 

dt  b 

F^-D  is  the  surface  flux  due  to  wave  breaking  minus  internal  dissipation  in 
the  mixed  layer.  E  is  the  rate  of  production  due  to  the  mean  shear  at  the 
base  of  the  mixed  layer.  TEM  models  neglect  E  and  attribute  the  time  rate  of 
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change  of  the  potential  energy  to  the  surface  flux  minus  internal 
dissipation.  DIM  models  do  the  opposite  and  consider  only  the  change  in 
potential  energy  due  to  mean  shear  at  the  base  of  the  mixed  layer. 

De  Szoeke  and  Rhines  (1976)  state  that  a  DIM  is  appropriate  for  the  first 
inertial  period  after  the  onset  of  a  wind  stress  forcing  a  fluid  initially  at 
rest.  For  longer  times  a  TEM  is  appropriate.  Klein  ana  Coantic  (1981) 
modeled  the  mixed  layer  using  a  second-order  turbulence  closure  model.  They 
concluded  that  the  significance  of  turbulent  kinetic  energy  input  due  to 
surface  wave  breaking  is  limited  to  the  first  few  meters  of  depth,  and  that 
when  the  depth  is  greater  than  this,  the  mixed-layer  deepening  is  driven 
primarily  by  the  shear  instability  at  the  base  of  the  mixed  layer.  This 
conclusion  supports  the  observations  of  mixed-layer  deepening  by  Price,  et  al . 
(1978)  who  found  that  the  deepening  is  highly  intermittent  on  the  storm  time 
scale  and  that  it  is  the  velocity  shear  at  the  base  of  the  mixed  layer  and  not 
the  magnitude  of  the  local  wind  stress  which  is  the  most  relevant  controlling 
factor. 


Because  of  the  short  time  scale  involved  related  to  a  storm  passage,  it 
is  clear  that  a  DIM  should  be  used.  A  TEM  could  only  be  valid  over  a  seasonal 
cycle,  although  its  validity  even  on  this  time  scale  is  not  certain  because 
mixed-layer  deepening  can  be  highly  intermittent  and  largely  confined  to  storm 
passages . 


The  first  slab  model  of  the  mixed  layer  using  the  DIM  mechanism  was  by 
Pollard,  et_  al_. ,  (1972).  It  is  based  upon  the  assumption  that  the  mixed  layer 
will  remain  marginally  dynamically  unstable.  This  assumption  is  manifested  by 
requiring  that  the  bulk  Richardson  number 


Ri  =  9Ap!U- 

p(u2+y2) 


(2) 


always  have  the  value  unity.  Here  h  is  the  mixed  layer  depth  and  Ap,  u,  and  v 
are  the  density  and  velocity  jumps  at  the  base  of  the  mixed  layer. 


As  an  alternative  to  requiring  that  Ri  remain  equal  to  one,  Phillips 
(1977)  postulated  that  the  rate  of  deepening  due  to  entrainment  is 


(u2+  v2)  F  (Ri) 


(3) 


where  F  is  an  unknown  function  of  Ri.  Price  (1979)  determined  F  using 
experimental  measurements  of  mixed  layer  deepening.  His  parametization  of 
entrainment  is  used  to  predict  mixed  layer  deepening. 

In  this  model  the  surface  front  moves  via  cross  stream  advection  in  the 
mixed  layer.  To  illustrate  this,  consider  initially  the  case  where  density  is 
conserved  in  the  mixed  layer.  Continuity  of  density  is  expressed  as 


dp  _  3p  , 

it  "3t 


v^P  =  0. 
3y 


(4) 


Here  v  is  the  cross  stream  (northward)  velocity  and  y  is  the  cross  stream 
(northward)  coordinate.  Equation  (4)  can  be  rearranged  to  the  form 
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f  \ 

9y 

3t 

v.  y 


P 


V. 


(5) 


This  form  of  (4)  shows  that  the  cross  stream  speed  holding  density  constant  is 
v.  Because  we  assumed  no  mixing,  the  surface  front  (or  any  surface  feature 
identifiable  by  its  density)  moves  with  speed  v. 


However,  in  general,  density  in  the  mixed  layer  is  not  conserved  due  to 
entrainment  of  cooler,  deeper  layer  water  at  the  base  of  the  mixed  layer.  The 
determination  of  the  density  balance  in  the  mixed  layer  and  the  speed  at  which 
the  surface  front  is  advected  requires  that  this  entrainment  be  described. 


II.  MODEL  DEVELOPMENT 

Our  model  consists  of  a  mixed  layer  with  a  narrow  density  front  overlying 
a  lower  layer  with  a  density  front.  Figure  1  illustrates  the  initial  state  of 
this  model  before  the  start  of  the  wind  stress  forcing.  We  start  with  the 
surface  front  overlying  the  subsurface  front.  The  density  field  varies  slowly 
on  each  side  of  the  surface  and  subsurface  fronts.  In  the  lower  layer  the 
downstream  velocity  U0  and  the  Brunt-Vaisiala  frequency  N  are  prescribed. 

They  vary  with  cross  stream  position  but  are  not  allowed  to  change  with 
time.  The  lower  layer  can  be  vertically  "pumped"  when  convergence  occurs  in 
the  mixed  layer.  At  the  base  of  the  mixed  layer  there  are  the  density  and 
velocity  jumps  Ap,  u,  and  v.  Initially  u  and  v  are  zero.  The  density 
jump  Ap  and  the  mixed  layer  depth  h  are  given  initial  values. 

A.  Mixed  Layer  Density  Balance 

When  entrainment  occurs, 

d£.=  9£+y3£_=^e_Ag.  (6) 

dt  3t  3y  3t  h 

where  he  is  the  increase  in  mixed  layer  depth  due  to  entrainment.  There  can 
also  be  a  flux  of  density  from  the  atmosphere  to  the  mixed  layer,  but  this 
effect  is  not  considered.  Price  (1981),  in  modelling  the  response  of  the 
upper  ocean  to  a  moving  hurricane,  found  that  air-sea  heat  exchange  plays  only 
a  minor  role  relative  to  entrainment  in  changing  the  mixed  layer  density. 

B.  Entrainment  Model 

In  order  to  predict  the  rate  of  mixed  layer  deepening  due  to  entrainment, 
we  use  a  DIM  (model)  by  Price  (1979).  Using  the  experimental  results  on 
stress-driven  entrainment  of  Kato  and  Philips  (1969)  and  Kanthe,  et  al . 

(1978),  Price  determined  that 

3he/3t 


5xlO"q'Ri"4  if  Ri  <c  1 
0  if  Ri  >  1 


(7) 
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C.  Mixed  Layer  Continuity 

The  equation  of  continuity  integrated  over  the  depth  of  the  mixed  layer 

^  +  -  wb  =  0  .  (8) 


Here  wb  is  the  vertical  velocity  (following  a  parcel)  at  the  base  of  the  mixed 
layer.  It  is  assumed  that  the  vertical  velocity  at  the  surface  is  zero. 

We  define  hc  to  be  the  change  in  depth  of  the  mixed  layer  due  to  a  cross 
stream  convergence  of  mass  in  the  mixed  layer.  Its  sign  is  positive  when  the 
increase  in  depth  is  positive.  Since 

<u °+u>i  (»> 

Equations  (8)  and  (9)  show  that 

3hc  ^  3((U„+u)h) 

at  ay  ax  u  ' 


The  cross  stream  convergence  must  be  balanced,  as  (10)  shows,  by  deepening  of 
the  mixed  layer  and/or  by  a  downstream  convergence.  Since  the  model  is  only 
two-dimensional  (2-D),  it  cannot  tell  how  the  cross  stream  convergence  is 
partitioned.  Therefore,  we  only  consider  the  two  limiting  cases  where 

3^  .  _  ^  .  »((Uq*0»>  . 

at  ay  ’  ax 

ahc  _  _  a((u0+u)h)  a ( vh )  . 

at  _  0  ’  ax  ay 

Using  the  same  terminology  as  Cushman-Roisin  (1982),  we  denote  (11)  to  be 
convergence  and  (12)  to  be  confluence. 


(ID 

(12) 


The  total  deepening  rate  ah/at  of  the  mixed  layer  is  the  sum  of  the 
deepening  rates  due  to  convergence  and  entrainment.  For  the  case  of 
confluence,  ah/at  is  equal  to  3he/at.  For  the  case  of  convergence 


3h  = 

at 


a(vh) 

'  ay 


if  Ri  £  1 


3(vh) 

ay 


if  Ri  >  1 


(13) 


D .  Jump  Condition  for  Density 

The  density  jump  across  the  interface  at  the  base  of  the  mixed  layer  is 

Ap  =  P£  -  P  (14) 
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where  p  is  the  mixed  layer  density  and  Pj£  is  the  density  in  the  lower  layer 
immediately  below  the  base  of  the  mixed  layer.  It  is  also  the  density  of  the 
water  entrained  into  the  mixed  layer. 

We  define  a  P£Q  as  the  initial  P&.  We  also  assume  that  the  lower  layer 
has  constant  stratification  N.  After  the  initiation  of  the  wind  stress, 
changes  as  the  mixed  layer  deepens.  Consider  the  case  of  confluence  where  the 
mixed  layer  deepens  only  due  to  entrainment.  For  this  case 


Ap  =  p£q  +  he  — —  -  p  (15) 

where  p0  is  a  reference  density  chosen  to  be  1  gm/cm^  for  simplicity. 

For  the  case  of  convergence  it  is  also  necessary  to  consider  the  pumping 
of  the  lower  layer.  The  initial  response  of  the  lower  layer  to  the  pumping  is 
solely  the  vertical  displacement  of  lower  layer  water  parcels.  This  simply 
moves  the  lid  of  the  lower  layer  up  or  down  leaving  p£  and  hence  the  validity 
of  (15)  unchanged.  Later  as  the  pumping  spins  up  the  lower  layer,  horizontal 
advective  terms  will  be  induced  there.  The  cross  stream  advection  of  density 
in  the  lower  layer  will  alter  p£  and  reduce  the  validity  of  (15).  However, 
because  we  are  modelling  the  mixed  layer  for  about  one  day,  we  ignore  this 
later  response  and  assume  that  (15)  remains  valid  for  the  case  of  convergence 
as  well. 


E.  Momentum  Balance 


Integrated  over  the  depth  h  of  the  mixed  layer,  the  momentum  balances  are 

-  (16) 


9(vh) 

9t 

where  we  have  assumed 


9h 


9t 


c  +  fuh  =  ^  +  Dy  -  &  f 


9P , 


p  \9y  9y 


f  u  =  -  1  !!° 

fU°  P  9y  » 

T(z=-h)  _  ^<5 

p  9t 


and 


9P  fo 
9y  "  9y 


qh9p  qAp  9h 
29y  "  P  9y 


(17) 


(18) 

(19) 

(20) 


In  (16)  and  (17)  D*.  and  Dw  are  terms  representing  frictional  damping.  Tx 
and  are  respectively  tne  east  and  north  components  of  the  wind  stress  T. 
Equation  (19)  used  by  Pollard,  et^  aj_.,  (1972)  says  that  the  stress  at  the  base 
of  the  mixed  layer  x(z=-h)  is  used  to  bring  the  momentum  of  entrained  lower 
layer  water  up  to  that  of  the  mixed  layer. 
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Outside  the  surface  front  the  frictional  damping  terms  are  probably  minor 
(especially  relative  to  the  wind  stress  t  which  we  assume  to  be  large). 

Inside  the  very  narrow  surface  front  they  may  not  be,  but  we  do  not  know  how 
to  accurately  represent  them.  Therefore,  two  different  types  of  simplifying 
assumptions  will  be  compared. 


The  first  simplifying  assumption  attempted  is  to  assume  Dx  and  Dy  to  be 
everywhere  negligible.  Equations  (16)  and  (17)  become 


3(uh) 

3t 


3 (vh)  _ 
3t 


3hc 

(f  - 

9Uo 

3u' 

j  vh 

3t 

V 

3y 

3y ' 

3hr 
c  + 

fuh 

»  i  - 

h/ 

'3P 

3t 

p 

P  \ 

,3y 

3Pr 


and 


(21) 

(22) 


At  the  start  of  each  run  of  the  model  before  the  advent  of  the  wind 
stress,  we  assume  that  the  mixed  layer  flow  is  in  geostrophic  balance,  i.e.. 


fu 


ghi  3p-j 

1  "  2p~  3y 


and 


(23) 


vi 


=  0- 


(24) 


The  i  subscripts  imply  initial  values.  The  mixed  layer  depth  h^  is  a 
constant  (3h1-/9y=0). 

The  alternative  type  of  simplification  of  (16)  and  (17)  stems  from 
considering  the  momentum  balances  away  from  the  surface  front.  Away  from 
surface  fronts  3p/3y-3p  /3y  is  a  minor  term  since  the  mixed  layer  density  and 
depth  gradients  are  small.  For  the  same  reason  u  should  be  much  less  than  U0 
so  that  3u/3y  can  be  neglected  relative  to  3lL/3y  in  (16).  As  before,  the 
frictional  damping  terms  are  neglected  as  well.  Hence,  away  from  surface 
fronts  (16)  and  (17)  become 


3(uh)  _  9hc 
3t  3t 


and 
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(26) 


The  first  type  of  simplification  of  (16)  and  (17)  assumes  that  the 
pressure  gradient  3P/3y-3P0/3y  in  the  surface  front  is  balanced  by  the 
inertial  and  Coriolis  terms.  However,  since  the  surface  front  can  be  a 
kilometer  or  less  in  width,  the  momentum  balance  there  is  likely  to  be  highly 
frictional.  Therefore,  the  second  type  of  simplification  of  (16)  and  (17) 
assumes  that  the  surface  front  is  in  frictional  balance,  i.e., 

3P  3Po  _  n 

37  '  3^r "  uy  (27) 
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This  assumption  is  very  convenient  because  equation  (26),  the  cross 
stream  momentum  balance,  is  now  appropriate  everywhere,  not  just  away  from 
surface  fronts.  Accordingly,  u  will  not  change  suddenly  from  inside  the 
surface  front  to  outside  it.  Therefore  it  is  consistent  to  assume  9u/9y  to  be 
minor  relative  to  9U0/9y  in  the  surface  front  as  well  and  also  assume  (26)  to 
be  valid  everywhere.  Using  (25)  and  (26)  to  describe  the  motion  of  the 
surface  front  makes  sense  if  one  thinks  of  the  surface  front  as  being  an 
arbitrarily  thin  interface  between  two  water  masses  of  differing  density.  The 
interface  moves  simply  with  the  speed  of  the  water  masses  immediately  on  each 
side  of  it. 

To  summarize,  two  different  simplifications  of  (16)  and  (17)  given 
respectively  by  (21)  and  (22)  and  by  (25)  and  (26)  are  tried.  The  major 
difference  between  the  two  is  that  the  first  type  of  simplification  includes 
the  pressure  gradient  term  3P/3y-9P0/3y  while  the  second  does  not.  For  later 
use  we  call  the  simpler  momentum  balance  which  excludes  the  pressure  gradient 
term  momentum  balance  A.  The  more  complete  momentum  balance  is  called 
momentum  balance  B. 

The  model  is  solved  numerically  using  a  100  meter  cross  stream  resolution 
(grid  scale)  and  a  100  second  time  step.  Momentum  balance  B  proved  to  be 
unstable.  To  counteract  this,  the  solution  was  smoothed  slightly  after  every 
time  step  with  the  intent  of  removing  grid  scale  instabilities.  This  was  done 
by  averaging  the  solution  in  each  grid  with  its  two  immediate  neighbors  using 
a  normalized  one-twelve-one  running  average.  The  model  solutions  assuming 
momentum  balance  A  and  confluence  were  obtained  with  and  without  smoothing. 
Except  for  slightly  broader  surface  fronts,  the  smoothing  made  little 
essential  difference.  For  better  comparison  with  the  solutions  assuming  B, 
all  solutions  assuming  A  displayed  in  the  results  are  those  that  have  been 
smoothed . 


III.  SIMULATING  HURRICANE  PASSAGE 

In  this  model  we  vary  the  wind  stress  in  direction  and  magnitude  with 
time  as  well  as  with  a  cross  stream  position  consistent  with  the  passage  of 
Hurricane  Dennis.  We  compute  the  wind  stress  by  assigning  the  storm's 

a.  wind  stress  distribution, 

b.  speed  and  direction,  and 

c.  initial  position  relative  to  the  2-D  section  of  the  surface  front. 

Because  the  size  and  strength  of  Dennis  were  only  roughly  known,  for 
simplicity  we  assume  the  storm  to  be  circularly  symmetric  and  to  have  a 
Gaussian  stress  envelope.  This  envelope  is  specified  by  assigning  the  wind 
stress  at  the  storm's  center  and  at  some  radius  from  its  center.  This  is 
analogous  to  the  manner  in  which  the  strength  of  a  hurricane  is  reported  by 
giving  the  wind  speed  at  its  center  and  the  radius  of  its  gale  force  winds. 
The  storm  may  have  any  starting  point,  speed,  and  direction,  but  the  speed, 
course,  and  strength  of  the  storm  is  held  fixed  during  its  passage.  During 
the  simulation  of  the  storm's  passage  the  direction  and  distance  to  the  storm 
as  a  function  of  cross  stream  position  are  continuously  computed.  These 
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parameters,  together  with  the  assigned  wind  stress  distribution  of  the  storm, 
allow  the  wind  stress  direction  and  magnitude  to  be  continuously  determined. 

Referring  to  figure  2  in  Horton  (1982),  we  see  that  Dennis  traveled  a 
path  essentially  parallel  to  that  of  the  Gulf  Stream  on  the  average  of  about 
40  km  to  50  km  to  the  south  of  the  subsurface  front.  The  figure  also  shows 
that  during  its  passage  Dennis  decreased  in  maximum  strength  and  was 
downgraded  to  a  Tropical  Storm.  However,  National  Hurricane  Center  advisories 
show  that  as  its  maximum  windspeeds  decreased,  its  radius  of  gale  force  (34 
knots)  winds  increased.  Based  upon  these  advisories  we  assume  the  wind  speed 
at  the  storm's  center  to  be  60  knots  (30  m/s)  and  the  radius  of  its  gaje  force 
winds  to  be  200  nm  (370  km).  Wind  stress  is  computed  from  wind  speed^using 

T  =  pa^D^10^  ^8) 


where  p,  is  air  density,  U^g  is  the  wind  speed  at  10m  above  the  sea  surface, 
and  Cn  is  the  drag  coefficient.  Bunker  (1976)  summarized  the  drag  coefficient 
calculations  of  many  investigators.  The  experimental  drag  coefficients  show 
much  scatter  and  are  evidently  not  known  accurately.  We  assume  Garrett's 
(1977)  form  for  Cp, 

CD  =  (0.73  +  0.069  U1Q)  x  10"3  (29) 

Using  this  expression,  Cp  equals  2.8x10“ 3  when  U^g  is  30  m/s,  and  1.9xl0"3 
when  U^g  is  17  m/s.  The  wind  stress  is  estimated  to  be  30  dynes/crrr  at  the 
storm's  center  and  to  be  6  dynes/crrr  at  a  radius  of  370  km. 

In  order  to  simulate  the  hurricane's  passage,  it  is  assumed  that  the 
storm  travels  due  east  44  km  to  the  south  of  the  initial  position  of  the 
surface  front  at  11  m/s  (40  km/hr).  We  also  try  two  other  storm  paths  to  show 
in  a  more  general  manner  how  a  hurricane  can  affect  the  Gulf  Stream  surface 
front . 


IV.  RESULTS 

The  three  different  storm  paths  used  to  model  surface  front  advection  are 
illustrated  in  figure  2.  For  each  of  these  paths  the  storm  travels  at  the 
same  speed,  11  m/s,  and  has  the  same  size  and  strength.  Figure  3  shows  the 
wind  stress  vectors  as  a  function  of  time  for  these  paths  at  the  initial 
location  of  the  surface  front  (the  location  of  the  subsurface  front).  In 
figure  4  the  displacements  of  the  initial  surface  front  are  shown  as  a 
function  of  time.  This  figure  is  computed  assuming  confluence  for  the 
continuity  balance  and  the  simpler  momentum  balance  A.  Essentially  the  same 
results  are  obtained  if  convergence  for  the  continuity  balance  or  the  more 
complete  momentum  balance  are  used  instead. 

Finally,  figures  5(a,  b,  c),  6(a,  b,  c),  and  7(a,  b,  c)  illustrate 
respectively  for  the  three  storm  paths  the  mixed  layer  densities  and  mixed 
layer  depths  as  a  function  of  cross  stream  position  at  several  times  (0,  8, 

15,  and  22  hours).  Figures  5a,  6a,  and  7a  assume  confluence  for  the 
continuity  balance  and  case  A  (where  the  pressure  gradient  3P/3y-9P0/3y  is 
neglected)  for  the  momentum  balance.  Figures  5b,  6b,  and  7b  also  assume 
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confluence  for  the  continuity  balance  but  in  conjunction  with  the  more 
complete  momentum  balance  B.  Figures  5c,  6c,  and  7c  use  the  simpler  momemtum 
balance  A  but  assume  convergence  for  the  mass  balance.  There  is  also  a  figure 
5d  which  is  identical  to  5c  except  that  the  prescribed  subsurface  front  in  6d 
is  much  weaker. 

A.  Storm  Path  One 


The  wind  stress  for  this  path  is  dominantly  to  the  west  as  shown  by 
figure  3.  Figure  4  shows  that  the  surface  front  is  advected  to  the  north  (to 
the  right  of  the  wind)  as  expected.  As  for  all  of  the  storm  paths,  there  is  a 
strong  inertial  oscillation  superimposed  on  the  net  displacement  of  the 
surface  front.  Figure  4  covers  half  an  inertial  cycle  in  which  we  see  the 
surface  front  oscillating  between  17.5  km  and  8.5  km  north  of  its  initial 
position.  We  assume  the  average  of  these  two  values,  13  km,  to  be  the  net 
displacement  of  the  surface  front  by  the  wind.  This  compares  with  a  shift  of 
18  km  observed  due  to  the  passage  of  Dennis  (Horton,  1982). 

Figures  5a,  5b,  and  5c  show  the  maximum  and  minimum  shifts  of  the 
surface  front  at  15  and  22  hours  respectively.  While  they  all  show  similar 
shifts  in  the  initial  surface  front,  this  similarity  does  not  extend  to  mixed 
layer  depth  and  density.  Figures  5a  and  5b  which  both  assume  confluence  give 
similar  mixed  layer  depths  and  densities.  However,  figure  5c  which  assumes 
convergence  differs  substantially  in  this  regard. 

Comparing  the  case  assuming  convergence  with  the  cases  assuming 
confluence,  the  mixed  layer  density  and  depth  at  8  hours  are  very  similar.  At 
15  hours,  though,  the  convergence  case  shows  strong  changes  in  mixed  layer 
depth  and  density  arising  from  cross  stream  convergences  of  mass.  From  about 
0  km  to  -10  km  there  is  upwelling  while  from  about  -10  km  to  -20  km  there  is 
downwelling.  The  upwelling  and  downwelling  are  due  to  wind  stress  induced 
convergences  of  mass. 

In  order  to  understand  why  this  occurred,  consider  equation  (25) 
where  the  local  time  rate  of  change  of  uh  is  neglected. 


Rearranging  this  equation  and  combining  it  with  (11)  yields 


_9_ 

3y 


f-3U0/8y 


(31) 


In  the  subsurface  front  3U0/3y  changes  much  more  rapidly  with  cross 
stream  position  than  Ty.  Thus 


^c  Tx _  d  UQ 

3t  "  (f-3U0/3y)2  3y2 


(32) 
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Equation  (32)  tells  us  that  a  wind  stress  with  no  curl  can  induce  intense 
Ekman  pumping  in  a  frontal  zone.  This  effect  was  first  described  by  Niiler 
(1969). 


Figure  3  shows  us  that  shortly  before  15  hours  Tx  was  large  and 
negative.  Between  0  km  and  -10  km  %d\}0/tyd  is  positive  while  between  -10  km 
and  -20  km  92UD/9y2  is  negative.  Therefore,  equation  (32)  predicts  upwelling 
(h  decreases)  between  0  km  and  -10  km  and  downwelling  between  -10  km  and 
-20  km.  This  is  exactly  as  figure  5c  shows  at  15  hours. 

After  16  hours  the  mixed  layer  water  parcels  inertially  oscillate 
with  little  wind  driving.  However,  adjacent  water  parcels  have  different 
inertial  amplitudes  and  frequencies  largely  because  of  cross  stream  changes  in 
the  effective  Coriolis  frequency  f-3LL/3y.  The  different  inertial  frequencies 
force  adjacent  fluid  parcels  to  have  their  respective  inertial  cycles  become 
increasingly  out  of  phase  with  time.  The  net  result  of  the  different  inertial 
amplitudes  and  phases  of  adjacent  fluid  parcels  is  intense  variability  at  the 
inertial  frequency  in  mixed  layer  mass  convergences  and  hence  mixed  layer 
depths . 


At  22  hours  the  differences  between  figures  5b  and  5c  have  become 
even  more  extreme.  Figure  5c  shows  that  by  this  time  a  dense  region  has 
formed  to  the  south  of  the  original  surface  front.  The  southern  edge  of  this 
dense  region  is  a  new  surface  front.  The  dense  region  has  formed  because  the 
large  gradients  in  mixed  layer  depths  due  to  assuming  convergence  have  led  to 
a  large  value  of  v9p/3y  there.  The  mixed  layer  depths  at  22  hours  in  figure 
5c  are  very  chaotic.  Since  after  16  hours  the  zonal  component  of  the  local 
wind  stress  is  weak,  the  chaotic  mixed  layer  depths  are  not  wind  induced  in 
the  manner  of  (32).  The  mixed  layer  depth  changes  are  similar,  though,  in 
that  they  are  for  the  most  part  due  to  large  cross  stream  changes  in  9U0/9y. 

In  order  to  verify  this,  another  run  of  the  model  assuming 
convergence  was  made  with  the  maximum  9lL/9y  reduced  from  -.6  (as  illustrated 
in  figure  3)  to  -.1.  The  result,  shown  in  figure  5d,  is  substantially 
equivalent  to  assuming  confluence. 

However,  the  results  assuming  convergence  are  at  least  superficially 
similar  to  observations  of  the  Iceland-Faeroe  Current  by  Von  Zweck  and  Tieg 
(1982).  They  observed  alternating  warm  and  cool  regions  (and  multiple  fronts) 
across  the  current.  In  this  area  the  horizontal  temperature-salinity  (TS) 
curves  were  the  same  as  the  vertical  TS  curves.  This  likeness  is  consistent 
with  the  alternating  warm  and  cool  regions  and  associated  fronts  being  due  to 
downwelling  or  upwelling. 

B.  Storm  Path  Two 


The  next  storm  path  is  also  parallel  to  the  stream,  but  is  north  of 
the  subsurface  front  the  same  distance  path  one  is  to  the  south.  This  path 
yields,  as  figure  3  shows,  meridional  wind  stresses  the  same  as  for  the  first 
storm  path,  but  zonal  wind  stresses  opposite  in  sign.  As  figure  4  shows,  the 
surface  front  initially  moves  to  the  north,  then  to  the  south.  The  net 
southern  displacement  is  expected  because  the  zonal  component  of  the  wind 
stress  is  positive.  Assuming  confluence  or  convergence  yields  essentially 
identical  surface  front  displacements. 
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As  with  the  first  storm  path,  the  two  cases  assuming  confluence  are 
very  similar.  At  15  hours  the  original  surface  front  is  near  4  km  and  has 
weakened  somewhat.  Near  -4  km  a  second  wider  surface  has  formed.  It  is  due 
to  entrainment  of  cooler,  underlying  water  to  the  south  of  the  subsurface 
front  than  to  the  north  of  the  subsurface  front. 

Figure  4  shows  that  after  13  hours  the  original  surface  front  is 
advected  to  the  south.  When  the  flow  is  to  the  south  in  the  mixed  layer, 
relatively  dense  mixed  layer  water  is  moved  over  suddenly  less  dense  lower 
layer  water.  This  leads  to  a  much  reduced  bulk  Ri  number  and  greatly  enhanced 
entrainment  at  the  base  of  the  mixed  layer.  Consequently  at  22  hours  we  see  a 
very  deep  mixed  layer  immediately  to  the  south  of  the  subsurface  front. 

At  22  hours  figure  4  shows  the  northern  edge  of  the  original  surface 
front  to  be  near  -18  km.  It  has  become  spread  out  and  diffuse  and  merges  into 
a  broad  decrease  in  mixed  layer  density  south  of  the  subsurface  front.  After 
22  hours  the  original  surface  front  loses  its  identity  and  becomes  impossible 
to  follow.  Figure  6c  which  assumes  convergence  also  shows  the  surface  front 
becoming  broad  and  diffuse  at  22  hours.  However,  an  exceedingly  deep  mixed 
layer  in  excess  of  200  meters  has  formed  over  the  subsurface  front.  As  with 
the  first  storm  path  this  is  primarily  due  to  sharp  cross  stream  changes  in 
3|J0/9y  there. 

C.  Storm  Path  Three 

For  this  final  situation  the  storm  path  is  purely  meridional  and  to 
the  east  of  the  section  of  the  modeled  current.  As  figure  3  shows,  this 
yields  wind  stresses  whose  meridional  components  are  always  negative.  Because 
the  winds  have  zonal  symmetry  during  the  storm's  passage,  one  might  expect  the 
surface  front  to  inertially  shake  about  its  original  position.  However, 
figure  4  shows  that  there  is  a  net  northward  displacement  of  the  surface  front 
which  occurs  because  the  mixed  layer  depth  increases  with  time.  This 
deepening  of  the  mixed  layer  dampens  the  movement  of  the  surface  front  making 
the  inertial  oscillation  nonsymmetric. 

Overall,  this  storm  path  gives  results  very  close  to  those  of  the 
first  storm  path.  The  major  difference  is  that  the  surface  front  at  22  hours 
in  figure  7a,  the  case  which  assumes  confluence  for  the  continuity  balance  and 
the  more  complete  momentum  balance,  has  lost  its  identity.  The  surface  front 
in  figure  7c  at  22  hours  has  weakened  substantially  as  well.  Its  remnant  is 
seen  at  5  km.  Of  interest  is  the  very  intense  new  surface  front  which  formed 
to  the  south  of  the  original  one.  Associated  with  this  new  surface  front  we 
again  see  large  subsurface  mixed  layer  depth  changes  due  to  cross  stream 
convergences  and  divergences  of  mass. 


V.  CONCLUSIONS 

Assuming  confluence  for  the  mixed  layer  continuity  balance,  two  degrees 
of  simplification  of  the  mixed  layer  momentum  balances  were  employed.  The 
more  complete  momentum  balance  includes  the  cross  stream  pressure  gradient  due 
to  cross  stream  changes  in  mixed  layer  density  and  depth  while  the  simpler 
momemtum  balance  does  not.  Including  or  disregarding  this  pressure  gradient 
made  little  difference  upon  frontal  displacement.  The  simpler  momentum 
balance  was  also  tried  assuming  convergence.  Again  the  observed  frontal 
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displacements  remained  essentially  the  same,  although  significant  differences 
in  the  details  of  the  mixed  layer  density  structure  and  depth  were  observed. 
This  observation  is  consistent  with  the  results  of  Price  (1981)  who  found  that 
mixed  layer  velocities  during  a  hurricane's  passage  are  primarily  sensitive  to 
the  wind  stress. 

The  first  storm  path  was  intended  to  roughly  simulate  the  passage  of 
Hurricane  Dennis  during  a  Gulf  Stream  survey  by  Horton  (1982).  The  model 
predicted  a  shift  of  the  surface  front  relative  to  the  subsurface  front  of. 
about  13  km.  This  compares  to  an  observed  shift  of  about  18  km.  Considering 
the  crudeness  of  our  knowledge  of  the  wind  stress  during  Dennis'  passage,  this 
agreement  is  quite  reasonable.  Also  obscuring  the  difference  between  the 
observed  and  predicted  shifts  is  the  inertial  oscillating  of  the  predicted 
surface  front  displacement. 

The  details  of  the  density  structure  in  the  mixed  layer  are  very 
sensitive  to  whether  confluence  and  convergence  are  assumed.  As  was  shown, 
these  greatly  different  results  arise  primarily  from  large  cross  stream 
changes  in  9U0/9y  in  the  subsurface  front.  The  changes  in  9U0/9y  induce 
substantial  cross  stream  convergences  in  mass  flux.  In  the  cases  assuming 
convergence  for  their  continuity  balances,  the  mass  convergences  cause  Ekman 
pumping.  These  cases  give  the  impression  of  being  too  extreme  and  thus 
unrealistic.  However,  they  strongly  suggest  that  substantial  variability  in 
mixed  layer  depths  at  the  inertial  frequency  should  be  expected  in  frontal 
regions . 

The  model  also  shows  that  the  rate  of  mixed  layer  deepening  due  to 
entrainment  is  strongly  dependent  upon  the  sign  of  the  cross  stream  advection 
in  the  mixed  layer.  Northward  advection,  as  associated  with  the  first  storm 
path,  results  in  relatively  light  mixed  layer  water  moving  over  suddenly  more 
dense  lower  layer  water  as  it  crosses  over  the  lower  layer  front.  A 
relatively  large  bulk  Richardson  number  and  a  relatively  low  entrainment  rate 
are  a  consequence.  Southward  advection,  as  occurring  with  the  second  storm 
path,  results  in  mixed  layer  water  moving  over  suddenly  less  dense  lower  layer 
water.  As  figures  6a  and  6b  show,  a  very  deep  mixed  layer  immediately  to  the 
south  of  the  subsurface  front  is  formed  as  a  consequence  of  the  low  bulk 
Richardson  number. 

The  primary  purpose  of  the  model  is  to  predict  the  displacement  of  an 
existing  surface  front.  The  model  is  judged  suitable  for  this  purpose.  Its 
most  pragmatic  version  is  the  version  using  the  simplest  momentum  balance  and 
assuming  confluence.  The  model,  however,  is  not  suitable  for  predicting  the 
details  of  the  mixed  layer  structure  when  a  strong  subsurface  front  is 
present.  A  two  dimensional  model  simply  cannot  distinguish  between  the  two 
limiting  cases  of  confluence  or  convergence  for  the  continuity  balance.  A 
three  dimensional  model  may  be  necessary  if  the  details  of  the  mixed  layer 
structure  in  a  frontal  zone  are  to  be  accurately  predicted. 
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CROSS  STREAM  POSITION  (km) 


Figure  1.  Initial  state  of  the  model  before  the  start  of  the  wind 

stress.  After  the  advent  of  the  wind  stress  the  mixed  layer 
density  and  the  density  at  the  top  of  the  lower  layer 

evolve.  The  lower  layer  stratification  N  and  the  lower  layer 
cross  stream  velocity  shear  9U0/9y  remain  fixed. 
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Figure  2.  The  three  storm  paths  used  relative  to  the  cross  section  of  the 
current  modeled.  Each  of  the  three  paths  are  straight  along 
which  the  storm  moves  with  a  constant  speed  of  11  m/s.  The 
first  storm  path  is  parallel  to  the  current  (eastward)  and 
44  km  to  the  south  of  the  initial  position  of  the  surface  front 
(the  position  of  the  subsurface  front).  The  second  path  is 
also  parallel  but  is  instead  44  km  to  the  north.  The  third 
storm  path  is  perpendicular  to  the  current  (northward)  and 
44  km  to  the  east  of  the  cross  section  modeled. 
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Figure  3.  The  wind  stress  vectors  versus  time  obtained  using  the  storm 
paths  described  in  figure  2  and  the  hurricane  model  described 
in  the  text.  These  wind  vectors  are  those  at  the  initial 
position  of  the  surface  front.  At  any  given  time  the  wind 
vectors  vary  with  cross  stream  position. 
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Figure  4.  Plots  of  surface  front  displacement  versus  time  since  start  of 
the  wind  stress  for  each  of  the  three  storm  paths.  The  wind 
stresses  used  are  illustrated  in  figure  2.  Positive  and 
negative  displacements  refer  to  north  and  south  displacements 
respectively. 
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Figure  5a.  Mixed  layer  density  and  depth  changes  due  to  storm  path  1. 

Assumes  confluence  for  the  continuity  balance  and  momentum 
balance  A. 
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CROSS  STREAM  POSITION  (km) 

Figure  5b.  Mixed  layer  density  and  depth  changes  due  to  storm  path  1. 

Assumes  confluence  for  the  continuity  balance  and  momentum 
balance  B.  (Same  as  figure  5a  except  assumes  momentum 
balance  B.) 
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CROSS  STREAM  POSITION  (km) 

Figure  5c.  Mixed  layer  density  and  depth  changes  due  to  storm  path  1. 

Assumes  convergence  for  the  continuity  balance  and  momentum 
balance  A.  (Same  as  figure  5a  except  assumes  convergence  for 
the  continuity  balance.) 
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CROSS  STREAM  POSITION  (km) 

Figure  5d.  Same  as  figure  5c  except  that  the  cross  stream  velocity  shear 
in  the  subsurface  front  is  much  weaker. 
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CROSS  STREAM  POSITION  (km) 

Figure  6a.  Mixed  layer  density  and  depth  changes  due  to  storm  path  2. 

Assumes  confluence  for  the  continuity  balance  and  momentum 
balance  A.  (Same  as  figure  5a  except  uses  storm  path  2  instead 
of  storm  path  1.) 
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Assumes  confluence  for  the  continuity  balance  and  momentum 
balance  B.  (Same  as  figure  5b  except  uses  storm  path  2  instead 
of  storm  path  1.) 
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CROSS  STREAM  POSITION  (km) 

Figure  6c.  Mixed  layer  density  and  depth  changes  due  to  storm  path  2. 

Assumes  convergence  for  the  continuity  balance  and  momentum 
balance  A.  (Same  as  figure  5c  except  uses  storm  path  2  instead 

of  storm  path  1.) 
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CROSS  STREAM  POSITION  (km) 

Figure  7a.  Mixed  layer  density  and  depth  changes  due  to  storm  path  3. 

Assumes  confluence  for  the  continuity  balance  and  momentum 
balance  A.  (Same  as  figure  5a  except  that  storm  path  3  is  used 
instead  of  storm  path  1.) 
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CROSS  STREAM  POSITION  (km) 

Figure  7b.  Mixed  layer  density  and  depth  changes  due  to  storm  path  3. 

Assumes  confluence  for  the  continuity  balance  and  momentum 
balance  B.  (Same  as  figure  5b  except  that  storm  path  3  is  used 
instead  of  storm  path  1.) 
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CROSS  STREAM  POSITION  (km) 

Figure  7c.  Mixed  layer  density  and  depth  changes  due  to  storm  path  3. 

Assumes  convergence  for  the  continuity  balance  and  momentum 
balance  A.  (Same  as  figure  5c  except  that  storm  path  3  is  used 
instead  of  storm  path  1.) 
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